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Abstract 

In this paper we investigate an integrable loop model and its connection with a su- 
persymmetric spin chain. The Bethe Ansatz solution allows us to study some properties 
of the ground state. When the loop fugacity q lies in the physical regime, we conjecture 
that the central charge is c = q — 1 for q integer < 2. Low- lying excitations are exam- 
ined, supporting a super diffusive behavior for q = 1. We argue that these systems are 
interesting examples of integrable lattice models realizing c < conformal field theories. 
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In Statistical Mechanics the basic difference between an ordinary local model ( vertex mod- 
els, spin chain systems ) and a loop model is that for the latter the total weight configuration 
cannot be written as a product of local weights. Being intrinsically non-local, loop models ap- 
pears as ideal paradigm for studying statistical properties of extended objects such as polymers 
( see e.g. refs. [|l| ). 

In this letter we investigate some critical properties of an integrable intersecting loop model 
in a two dimensional square lattice 0. The fact that intersections between the polygons 
configurations are allowed makes this loop model very interesting. In this case it is not clear at 
all how to find a transformation to a standard local model 0, ||. The model can be considered 
as a Lorentz gas of particles moving through a set of scatters randomly distributed in the 
nodes of the two dimensional square lattice [|, |[ 0]. The scatters are double-sided mirrors 
allowing right-angle reflections and they are placed along the diagonals of the square lattice. 
The particles move along the edges of the lattice and arriving on a node can be scattered to 
the left, to the right or pass freely in the case of absence of a scatter. We denote by w a , Wb and 
w c the Boltzmann weights corresponding to these three possibilities, respectively. By imposing 
periodic boundary conditions, each particle follows a closed path. If for every closed loop we 
assign a fugacity q, then the partition function Z can be written as 



where n a ,n& and n c are the number of weights w a , Wb and w c appearing in a configuration, 
respectively. We noticed that strictly when w c = the closed loops configurations no longer 
intersect. In this limit the partition function (1) can be seen as a graphical representation of 



Despite the inherent non-locality of this loop model it is still possible to formulate a purely 
local condition that two different transfer-matrices commute for arbitrary system size. This is 
a sufficient condition for integrability and it imposes a restriction on the manifold of possible 
weights w a , Wb and w c configurations. It turns out that the intersecting loop model is integrable 
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the critical g 2 -state Potts model 0]. 
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|2| if the Boltzmann weights are parametrized as follows 

q 

w a (fj) = 1 - /I, W b (fj) = fi, w c (ii) = (1 - -)/i(l - //) (2) 

Clearly, the interesting physical regime is when q < 2, where all the "probabilities" w a , Wb 
and Wb can be made positive for < \i < 1. Anyway, it is interesting to noticed that when g 
is an integer > 2 these weights reproduce ( after a rescaling of fi ) the ones appearing in the 
S- matrix solution of the 0(q) invariant non-linear sigma-model ||. This strongly suggest that 
a mapping onto a ordinary local vertex model should not be completely ruled out, at least for 
some values of q. Indeed, here we argue that when q G Z the intersecting loop model can be 
realized in terms of a local spin chain which is invariant by the superorthosympletic Osp(n\2m) 
superalgebra. This is an important observation, since it will allow us to study the physical 
regime of (2) for q integer < 2. We then use the fact that the Osp(n\2m) super spin chain 
is solvable by the Bethe Ansatz in order to find that their critical properties are governed 
by c < conformal field theories. We remark that c < theories appear to be useful in 
condensed matter systems such as the quantum Hall effect || |10| , disordered models and 



polymer field theories ||12|| . This means that our models could be relevant lattice paradigms for 



some condensed matter applications, and being integrable, they might provide further physical 
insight as well. 

Essential to our approach is to observe that the weights (2) can be derived in the context 
of a standard Yang-Baxter solution for a local vertex model. These weights are in one-to-one 
correspondence to the generators of a degenerated point of the Birman-Wenzel-Murakami 
algebra ||13|| . This algebra is generated by the identity Ii, a braid hi and a Temperely-Lieb 
operator Ei acting on sites i and % + 1 of a quantum spin chain of length L. On the degenerate 
point the braid operator becomes a generator of the symmetric group, namely 

hk±xk = b i ± 1 b i b i±1 , b\ = I h bibj = bjbi if \i - j\ > 2 (3) 

and the other set of relations closing the degenerated point of the braid-monoid algebra ( see 
e.g ref. |0J ) are 

EiE i±1 Ei = E h El = qE u E l E j = E^ if \i - j\ > 2 (4) 



and 

biEi = Eibi = Ei, Eibi±\bi = bi±±biEi±i = E{Ei±\ (5) 



It is not difficult to see that relations (3-5) can be baxterized [|14], [15| , providing us with a 
rational solution of the Yang-Baxter equation having precisely the weights w a , Wb and w c . To 
make further progress we search for a representation of the algebraic relations (3-5). At least 
for integer q, such representation can be found in terms of the invariants of the superalgebra 
Osp(n\2m) [|14| , 15| . This superalgebra combines the 0(n) symmetry and the simpletic Sp(2m) 



algebra, and the integers n and 2m play the role of the number of bosonic and fermionic degrees 
of freedom ( see e.g. ref. [|l6j ). The braid operator bi becomes the graded permutation between 



the (n + 2m) degrees of freedom which is defined by [17 

n+2m 



h= {-^Y {a)p{h) e a c®e hd (6) 

a,b=l 



where p(a) is the parity distinguishing the bosonic (p(a) = for a — 1, • • • , n) and the fermionic 
(p{a) = 1 for a — n + 1, ■ ■ ■ ,n + 2m) elements. Explicit expressions for the monoids Ei have 
been recently discussed in detail in ref. |15| . The important point which matters here, however, 
is that the fugacity q is precisely the difference between the number of bosonic and fermionic 
degrees of freedom. More precisely, we have 

q = n — 2m (7) 

The formulation of the corresponding transfer-matrix has to respect the bosonic and the 
fermionic graduations |17|, [18| . This is possible by writing T(X) as the supertrace of an auxiliary 
monodromy operator, T(A) = ^2(—l) p ^T aa , where A stands for the horizontal space of 
(n + 2m) variables of the vertex model. As usual the monodromy matrix is composed by a 
product of vertex operators C^iW which are given by 

C Ai (X) = (1 - I - \)h + A(l - I - \)h + XEi (8) 

Performing the scale A — > — q/2) in equation (8), it is straightforward to see the 
correspondence between the weights (2) and the operators Jj, bi and E^ The corresponding 



local Osp(n\2m) spin chain Ti is proportional to the logarithmic derivative of the transfer 
matrix around the regular point A = where the vertex operator reduces to the graded 



permutation |I7|, [0| . The Hamiltonian can then be expressed in terms of the operators bi and 
Ei as 

_l r i i 

(9) 



where the sign in (9) is chosen to select the antiferromagnetic regime of the theory. This super- 
symmetric Hamiltonian, with periodic boundary conditions imposed, admits a Bethe Ansatz 
solution . This means that the eigenenergies E(L) on a ring of size L are parametrized in terms 
of complex set of variables {A"}, satisfying coupled non-linear Bethe Ansatz equations. These 
equations are equivalent to the analyticity of the eigenvalues of the corresponding transfer 
matrix and also reflect the underlying Osp(n\2m) group symmetry. By taking into account 
these properties ( analytical Bethe Ansatz approach) we can derive that they are given by 
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and the eigenenergies are parametrized by {A!- } 
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E(L) = -J2—n\ + L (ii) 

fe[AS x) ] 2 + l/4 

where C a b is the Cartan matrix, r is the number of roots and rj a is the normalized length 
of the a-th root of the Osp(n\2m) superalgebra. We recall that the Dynkin diagrams of the 
models with n — 1 are special. Consequently a peculiar two-body phase shift occurs for the 
last root {A| m }. For an algebraic Bethe Ansatz derivation of equations (10,11) for some classes 
of Osp(n\2m) models as well as for further technical details we refer to ref. ||15|| . 

We now turn to the study of the critical behaviour of the super spin chain (9). The 
existence of a Bethe Ansatz solution allows us, in principle, to calculate the eigenvalues E(L) 
for quite large values of L, providing us with reasonable estimates of the finite size effects. For 
a conformally invariant system, the universality can then be determined by exploiting a set 
of important relations satisfied by the eigenvalues on a strip of size L JT^]. For example, the 



central charge c is related to the ground state energy E (L) by |2 

E (L) 



L 00 6L 2 



+ 0(L- 2 ) (12) 



where eoo is the ground state energy per particle in the thermodynamic limit and v s is the 
sound velocity These two parameters can be determined exactly from the unitarity and the 
crossing properties ( around A = l — q/2 ) of the transfer matrix T(A). In fact, these properties 
together imply that, in the thermodynamic limit, the largest eigenvalue A (A) of the transfer 
matrix satisfies the relations 

A (A)A (-A) = [1 - A 2 ][(l - q/2) 2 - A 2 ], A (A) = A (l - q/2 - A) (13) 

Solving this equations with the restriction that the solution is free of zeroes in the physical 
strip 0<A<1 — q/2 and taking its logarithmic derivative at A = 0, we find that 

e - = -1^72 H + 2^> - %b> + 21n(2) ) + 1 (14) 

where ip(x) is the Euler function. The sound velocity measures how the energy scales with 
the low momenta. If we recall that equations (13) are identical to the one we solve to find the 
crossing factors in a relativistic S"-matrix theory, we can obtain that the appropriate relativistic 
scale is given by 

«■ = T^f2 < 15 > 

We now have the basic ingredients to begin our analysis of the finite size effects for the 
ground state energy. Let us first consider the case when the fugacity is one . For this value 
of q, the partition function of the intersecting loop model is trivial (constant) and therefore 
Eq(L)/L is exactly for any size L. In other words, all the finite size corrections to the 
ground state are null, and in particular c = 0. However, from the spin chain point of view, 
this scenario is far from being trivial, and provides us with an important check concerning the 
loop model «-> spin chain mapping. The simplest spin chain giving us q = 1 is the Osp(3\2) 
model. Its spectrum is given in terms of one level nested Bethe Ansatz and the Bethe equations 
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are parametrized by two sets of variables {X^\ Xj}. The ground state is characterized by a 
complex root distribution, forming "fractional" strings of the following type 

\V =$ 1) ±i/4 + 0(e- L ), Af=£f ±i/4 + 0( e - z ) (16) 

By solving numerically the corresponding Bethe Ansatz equation for some values of L and 
substituting the value of {A^ 1 } in equation (11), it is remarkable to see how the Bethe Ansatz 
roots conspire together in order to produce the simple result E (L) = — 3L exactly. Note that 
for this model that e ro = —3 (see equation (14)). Although similar effect has been observed 
before in fine tuned anisotropic spin chains ||21|| , to the best of our knowledge, this is the 
first time that such simplification is noted in a free-parameter (isotropic) set of Bethe Ansatz 
equations. This gives confidence to investigate the super spin chains for other values of q < 2. 

For q 1, the procedure described above can also be used, and we have analysed the 
equations (10,11) for sizes up to L = 80. In table 1 we show our estimates for the central charge 
cfor the Osp(2\2), Osp(l\2), Osp(l|4) and Osp(l\6) supersymmetric spin chains, corresponding 
to the values q = 0, —1, —3, and —5, respectively. In our numerical analysis we already 
have considered the presence of logarithmic contributions of 0(1 /L 2 ln(L)) to the finite size 
corrections of the ground state. We remark that this kind of correction typically cannot be 
overcome by standard transfer matrix or Hamiltonian diagonalization due to size limitations. 
All the results lead us to the following conjecture for the central charge behaviour for these 
models when q is a integer < 2 

c = q-l (17) 

This formula also reproduces the central charge in the limit q = 2. For this point, the 
weight w c vanishes and hence we expect the critical behaviour of the isotropic 6-vertex model. 
Furthermore, the ground state of the Osp(l\2n) models (q = 1 — 2n) are parametrized by 
real roots, and by using an analytical method developed in ref. |23| we obtain c = — 2n, in 
accordance with equation (17). Finally, it is interesting to note that formula (17) can be derived 
in the context of a super Sugawara construction of the Virasoro algebra developed by Goddard 



et al [24J]. Following this work the central charge at certain level k is c = ksdimG / (k + Qg/2), 



where sdimG and Qg is the superdimension and the value of the quadratic Casimir of the 
superalgebra G. This latter expression applied to the Osp(n\2m) superalgebra with k — 1 
gives c = n — 2m — 1 , which agrees with the formula we proposed for the central charge of the 
loop model when q = n — 2m < 2. 

In principle we can use a similar approach to study the excitations via numerical or an- 
alytical analysis of the Bethe Ansatz equations (10,11). This study is of particular physical 
relevance for the Lorentz lattice gas q = 1 model . In this case an interesting quantity is the 
fractal dimension df of the loops, which characterizes the diffusion properties of the particles 
through the lattice |7|. Recent numerical simulations |F| predicted a superdiffusive behaviour 
df = 2 as long as the density of mirrors is smaller than one (w c ^ 0). To lend some theoretical 
support to this class of universality we have studied the finite size corrections for the lowest 
excitation present in the q = 1 model. We find that this excitation is of spin wave type, made 
by taking out one root A| from the ground state configuration. In table 2 we present the finite 
size sequences for the exponent x = 2h, where h is the conformal weight. We see that the 
extrapolated value for h is very small, indicating the presence of a zero conformal weight and 
consequently predicting df = 2 — 2h = 2 whithin reasonable precision. We remark that better 
numerical data was prevented by strong logarithmic corrections []7| . 

A second example where the study of excitations plays an important role is for the su- 
persymmetric Osp(l\2) model (g = — 1). Its spectrum is parametrized by a single non-linear 
equation and some low-lying states are amenable of analytical study (real Bethe Ansatz roots). 
This model has c = —2 and to obtain the Ramond sector, having lowest conformal dimension 
h = —1/8, we need to impose antiperiodic boundary condition on the even U(l) sectors of 
the superspin chain ||22|| . The twisted Hamiltonian still has a continuous SU(2) symmetry and 
the U(l) charge plays the role of a fermionic index . In principle, there is no inconsistency to 
interpret the twisted sector as an exited state, even though its energy lies below the ground 
state, because the super spin chain is indeed non hermitian (though the eigenvalues are real). 
This scenario is in remarkably similar to the one proposed in refs. [TI| to explain certain 
properties of edge excitations in the fractional Hall effect. Furthermore, concerning the eigen- 
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spectrum of the primary fields, the twisted spin chain could either be seen as a c = —2 theory 
or as c = 1 Coulomb gas with radius of compactification p = Vz. We also noticed that the 
Osp(2\2) spin chain is the prototype for realizing both c = — 1 and c = 2 field theories. Again, 
our numerical data supportes the presence of a field with dimension h = —1/8 in the twisted 
spin chain. Very recently, such relations between c < and c > have been proposed and 



discussed in the literature |25[. This suggests that the lattice models discussed in this paper 
are realizations of non-unitary (c < 0) conformal theories that are related ( via appropriate 
twisting ) to c > systems. In particular our numerical result for the twisted q = 1 model 
indicates central charge c = 3, suggesting that the continuum limit of the Lorentz gas might 



be governed by a N = 2 supersymmetric field theory ( see e.g. ref.[26 ) . It remains to be seen, 
however, what kind of physical information can be obtained from this observation. 

In summary, a solvable loop model has been mapped onto a super spin chain and we have 
found its central charge behaviour for integer values of the fugacity in the physical regime. We 
hope that our discussions lend support to the idea that these models appears to be the ideal 
lattice candidates for describing conformal properties of relevant condensed matter systems. 
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64 


-1.07179 


-1.95812 


-3.98721 


-6.04643 


72 


-1.06911 


-1.95973 


-3.98593 


-6.03771 


80 


-1.06683 


-1.96107 


-3.98508 


-6.03125 


Extr. 


-1.01 (1) 


-1.996 (2) 


-3.985 (3) 


-5.997 (1) 



Table 1 : Finite size sequences for the extrapolation of the central charge for the models 
Osp(2\2) (q = 0), Osp(l\2) (q = -1), Osp(l|4) (q = -3) and Osp(l\6) (q = -5). 

Table 2 : Finite size sequences for the extrapolation of the lowest conformal dimension 
x = 2h of the q — 1 model. 
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x = 2h 


16 


8.464651 10~ 2 


32 


6.764093 10- 2 


48 


6.032298 10~ 2 


64 


5.592822 10~ 2 


80 


5.288533 10~ 2 


96 


5.060282 10~ 2 


112 


4.88002 10~ 2 


128 


4.73245 10~ 2 


Extr. 


1.33 10~ 2 
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